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ABSTRACT 


I 


This  report  describes  three  phases  of  work  towards  improving  and 
|  extending  the  three-dimensional  ship  boundary  layer  calculation  method  of 

von  Kerczek  (1982)  developed  under  prior  GHR  support.  The  three  items  dis¬ 
cussed  in  this  report  are  the  following:  (1)  the  method  has  been  extended 
to  deal  with  roughness  of  the  ship  hull;  (2)  a  viscous  inviscid  interac- 
|  tion  scheme  has  been  developed  but  not  yet  computationally  implemented;  and 

(3)  the  computer  implementation  of  the  basic  von  Kerczek  momentum  integral 
method  has  been  improved  in  several  details.  These  are  described  in  this 
report. 


1.0 


INTRODUCTION 


Under  prior  GHR  support  von  Kerczek  (1982)  (henceforth  referred  to 
as  K)  developed  a  momentum  integral  ship  boundary  layer  calculation  method 
which  seems  to  be  a  promising,  relatively  simple  and  useful  way  of  estimating 
the  viscous  flow  and  some  of  its  consequences.  This  method  is  capable  of 
dealing  with  reversed  cross-flow  boundary  layer  velocity  profiles  and  thus 
can  be  made  more  acucrate  than  prior  momentum  integral  methods  for  ships 
(see  von  Kerczek  and  Langan,  1979). 

The  method  of  K  has  already  been  put  to  use  to  calculate  the 
total  drag  of  a  destroyer  hull  (von  Kerczek,  et  al . ,  1983)  and  as  an  aid  in 
determining  appendage  drag  (Stern  and  von  Kerczek,  1983).  These  applica¬ 
tions  are  still  tentative  and  require  further  verification  and  refinement. 
However,  the  results  are  encouraging  and  also  have  pointed  the  way  towards 
some  of  the  more  important  weaknesses  of  the  method  of  K. 

Under  further  GHR  support,  for  which  this  report  is  presented, 
three  items  have  been  considered  to  improve  the  K  method.  These  are  the 
fol lowing: 

(1)  The  method  of  K  has  been  extended  to  deal  with  roughness  of  the 
ship  hull . 

(2)  A  method  for  viscous-inviscid  interaction  has  been  developed, 
but  not  yet  computationally  implemented. 

(3)  The  computer  implementation  of  the  basic  method  of  K  has  been 
improved  in  several  details  making  it  more  reliable  and  accurate. 

This  report  details  each  of  these  efforts. 
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2.0  A  GENERAL  PROCEDURE  FOR  INCLUDING  ROUGHNESS  EFFECTS  IN  THE 

INTEGRAL  BOUNDARY  LAYER  EQUATIONS 

2.1  Introduction 

Most  roughness  analyses  are  based  on  Nikuradse's  sand-grain 
experiments  and  fitting  of  the  law-of-the-wall  velocity  profiles  to  experi¬ 
mental  data.  Several  correlations  have  been  proposed  to  relate  real 
roughness  heights,  spacings  and  geometries  to  an  equivalent  sand-grain 
height  so  that  Nikuradse's  data  can  be  used.  Examples  of  such  correlations 
can  be  found  in  Betterman  (1966),  Simpson  (1973)  and  Dirling  (1973).  These 
correlations  are  not  very  accurate  for  closely  spaced  roughness  elements 
such  as  those  found  on  a  ship's  surface.  The  correlations  in  these  refer¬ 
ences  are  biased  toward  data  for  roughness  elements  with  large  spacings. 

More  meaningful  methods  that  account  for  actual  roughness  effects 
are  those  employed  by  Christoph  and  Pletcher  (1983),  Finson  and  Clarke  (1980) 
and  Lin  and  Bywater  (1982).  These  techniques  calculate  the  form-drag  con¬ 
tributions  of  individual  elements.  Roughness  elements  are  assumed  to  occupy 
no  physical  space.  The  governing  boundary  layer  equations  are  cast  in  a 
form  to  account  for  the  blockage  effects  of  the  roughness  elements.  Terms 
that  act  in  the  streamwise  direction  are  multiplied  by  the  factor  [l-D(y)/£], 
where  D(y)  is  the  element  width  at  height  y  and  l  is  the  average  center-to- 
center  spacing.  Terms  that  act  in  the  direction  normal  to  the  streamwise 
direction  are  multiplied  by  [l  -  iTD2(y)/(4£2 )] .  The  effect  of  roughness  in 
the  momentum  equations  is  described  by  a  sink  term  and  represents  form  drag. 

To  date,  the  form  drag  approach  has  been  taken  only  with  a  finite- 
difference  solution  of  the  boundary  layer  equations.  In  order  to  apply  this 
method  directly  to  a  3-dimensional  integral  technique  such  as  that  of  K,  a 
skin  friction  correlation  that  explicitly  accounts  for  roughness  height, 
spacing  and  geometry  has  to  be  specified.  This  is  not  yet  possible;  instead, 
the  following  approach  is  taken  here. 


2-1 


Based  on  results  of  numerous  computer  runs  with  a  variety  of  rough¬ 
ness  shapes  and  spacings,  Finson  (1982)  suggested  a  correlation  relating 
equivalent  sand-grain  height  to  roughness  height,  spacing  and  geometry. 

This  correlation  was  shown  by  Finson  to  be  accurate  for  closely  packed  ele¬ 
ments.  Using  calculated  equivalent  sand-grain  heights,  a  method  was  devel¬ 
oped  in  this  study  to  calculate  skin  friction,  boundary  layer  thickness  and 
velocity  profiles  based  on  classical  sand-grain  roughness  relations.  The 
details  are  described  below. 


Solution  Procedure  for  2-D  Boundary  Layers 


It  is  desirable  to  include  roughness  in  the  full  3-D  boundary 
layer  equations  using  Head's  entrainment  technique  as  implemented  by  K.  To 
begin  the  development,  the  roughness  methodology  is  most  easily  discussed 
in  the  framework  of  a  2-D  flow.  Extension  to  3-D  and  its  implementation 
are  discussed  later.  The  2-D  momentum  integral  equation  is 


(2-1) 

(2-2) 

(2-3) 


The  skin  friction  coefficient  cf.  Head's  shape  factor  6,  and  the  rate  of 
entrainment  function  F(G)  are  modelled  in  the  same  way  as  in  K.  In  the 
present  study  it  is  assumed  that  roughness  effects  are  accounted  for  simply 
by  including  roughness  in  the  skin  friction  model.  The  momentum  thickness 
0  and  the  shape  factor  H  are  affected  by  roughness  through  integration  of 
equations  (2-1)  and  (2-2).  The  process  of  entrainment  is  mainly  controlled 
by  the  velocity  defect  in  the  outer  part  of  the  boundary  layer  and  therefore 
the  entrainment,  as  used  in  K,  should  not  be  affected  by  roughness.  The 
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magnitude  of  the  entrainment  variables  (G  and  F)  will  change,  however,  due 
to  their  dependence  on  H. 


Roughness  effects  are  modelled  in  c^  through  the  classical  velocity 
shift,  that  is. 


u*  =  u*  -  AB  (2-4) 

where 

u!  =  7  in  y+  +  CONSTANT,  K  =  0.41 

S  l\ 

A6  =  i  in  -  3 
and 

y+  ^  yu*/v,  =  kes  u*/v 
U+  =  u/u*,  u*  =  (tw/p)1/2. 

The  subscripts  s  and  r  refer  to  smooth  and  rough  conditions,  respectively. 
Since 


(2-5) 


(2-6) 


“edge  *  '  <2‘7> 

the  rough  wall  skin  friction  coefficient  Cfr  is  related  to  the  smooth  wall 
skin  friction  coefficient  c^s  by 

(2/cfr)1/2  =  (2/cf$)1/Z  -  AB.  (2-8) 

If  cfs  and  kgs  are  known,  then  cfr  is  easily  solved  from  equation  (2-8).  This 
is  done  by  rewriting  equation  (2-8)  as 

(2/c  r)1/2  -  £  In  (2/cfr)1/2  =  (2/cf$)1/2  -  i  £n  Rek  •»  3 
{Rekp„  =  Ue  kes/v) 


(2-9) 


and  using  an  iteration  scheme,  e.g.,  Newton-Raphson,  to  solve  for  cfr. 
method  used  for  determining  kes  is  discussed  next. 


The 


In  the  form  drag  analysis  outlined  in  the  introduction,  the  wall 
shear  is  given  by 


where 


f(y)  = 


1  -  D(y)A 
1  -  TrD2(y)/(4£2 ) 


(2-10) 

(2-11) 


and  Cq  is  a  form  drag  coefficient.  Taking  u  =  Up  =  constant,  cD  =  constant, 
and  evaluating  f(y)  at  y  =  k/2,  equation  (2-10)  becomes 


cfr  =  cfs  +  ^  cD  f(k/2)/Xk. 


(2-12) 


The  parameter 


-i 

Xk  =  projected  frontal  area  of  the  elements  per  (2-13) 

unit  surface  area. 


Based  on  numerous  calculations  for  c^  =  0.6,  Cfs  at  Re  =  10  ,  and  k/6  =  1, 
Finson  (1982)  recommends  the  correlation 


~  =  0.247  +  0.234  log  X. 

Ug  10  K 

So  for  the  above  conditions  equation  (2-12)  becomes 


(2-14) 


=  Ci 


(Re0  =  10s ) 


+  0.6(0.247+  0.234 


log  X.) 

10  K 


f(k/2)/Xk.  (2-15) 


A  simple  relation  can  be  obtained  for  the  equivalent  sand-grain  height,  ke$, 
if  Cf$  is  approximated  by 
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(2/Cf$)/z  =  in  ReQ  +  C  (flow  geometry,  pressure  gradient).  (2-16) 


Equation  (2-9)  can  be  rewritten  as 

(2/cfs)^  -  |  £n(2/cfr)  =  ^  £n(kes/0)  +  C  +  3. 


(2-17) 


A  general  expression  for  ke$  is  obtained  by  using  Finson's  correlation  (0=k) 
and  rearranging  equation  (2-17)  to 


=  exp  [K(C+3-  (2/c  *r)Ss  +  [  n(2/c*r)h]. 


(2-18) 


For  C  equal  to  its  flat  plate  value  of  5.5,  c  +  3  =  8.5.  This  represents 
flat  plate  flow  in  a  fully  rough  regime. 

2.3  Extension  to  3-D 

The  approach  taken  for  3-D  flows  is  a  direct  extension  of  the  2-D 
roughness  analysis  just  discussed.  The  streamwise  skin  friction  coefficient, 
cf ^ ,  that  occurs  in  the  streamwise  momentum  equation  is  replaced  by  the  rough 
wall  skin  friction  model  outlined  above.  Rough  wall  values  of  0^,  H  and  0£1 
(see  K  for  notation)  are  then  generated  by  solving  equations  (2-6)  -  (2-8)  of 
K.  The  cross  flow  skin  friction  coefficient,  c.^,  is  computed  from 


Cf2  =  Cfj  tan8 


(2-19) 


where  tan8  is  calculated  from  equation  (5)  of  K.  It  should  be  noted  that 
the  rough  wall  skin  friction  calculation  by  equation  (2-9)  requires  knowl¬ 
edge  of  the  smooth  wall  skin  friction.  Therefore,  the  smooth  and  rough 
wall  flows  need  to  be  calculated  in  parallel. 


In  order  to  obtain  the  crossflow  boundary  layer  thickness,  02j 
and  012»  it  is  necessary  to  know  how  roughness  affects  the  streamwise  and 
crossflow  velocity  profiles.  Roughness  effects  were  included  in  the  stream- 
wise  velocity  profile  as  follows.  Equation  (2-4)  is  written  as 


(2-20) 


ur/ur  =  us/us*  "  1  £n  ^kes  ur /v^  +  3 
This  is  rearranged  as 

Ur/Ue  =  Us/Ue  {cfr/cfs)H  +  <cfr/2^  C*  **((Z/cfr)'i)  (2-21) 

-  |  (ke$  ue/v)  +  3] 

where  us/Ug  is  a  smooth  wall  boundary  layer  velocity  profile.  For  the  results 
of  this  study,  the  profile  of  Kool  (see  K)  was  used.  The  crossflow  profile 
is  then  calculated  by  using  equation  (2-21)  as  the  streamwise  profile  in 
equation  (3)  of  K.  It  is  also  necessary  to  know  the  rough  wall  nominal 
boundary  layer  thickness.  It  is  calculated  from  the  relation 

6r  =  8n  (G  +  H)  (?-22) 

where  e^,  G  and  H  are  rough  wall  values. 

2.4  Computational  Results 

The  roughness  approach  described  above  was  validated  by  comparison 
to  the  2-D  flow  roughness  data  of  Schlichting  (1937)  and  Healzer,  Moffat  and 
Kays  (1974).  First,  comparisons  were  made  to  Schlichting's  channel  flow  data 
where  various  arrangements  of  roughness  elements  were  used.  Figure  2-1  gives 
skin  friction  coefficient  data  versus  distance  Reynolds  number  (as  presented 
by  Finson  and  Clarke  (1980))  for  spherical  segments  with  element  spacings 
of  L/D(o)  =  1,  2.5  and  5.  The  present  predictions  are  in  excellent  agree¬ 
ment  with  the  data  except  at  low  Reynolds  numbers  and  correctly  predict  the 
effect  of  element  spacing.  Figure  2-2  compares  the  rough  wall  velocity  pro¬ 
file,  equation  (2-21),  with  profile  data  taken  by  Schlichting.  the  compari¬ 
son  of  the  profiles  is  for  the  top  of  the  roughness  elements  out  towards 
the  edge  of  the  boundary  layer.  It  appears  that  equation  (2-21)  overpre¬ 
dicts  the  momentum  deficit  close  to  the  wall.  This  overprediction  is  due 
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to  the  presence  of  the  roughness  elements. 


The  incompressible  flow  experiments  of  Healzer,  et  al .  (1974)  pro¬ 
vide  data  on  skin  friction,  heat  transfer  and  velocity  profiles  for  rough 
surfaces  with  and  without  blowing.  The  rough  surface  consisted  of  a  regular 
array  of  tightly  packed  (L/D(o)  =  1)  hemispherical  elements,  each  0.05  inches 
in  diameter.  Tests  were  run  at  freestream  velocities  of  32,  139  and  242 
ft/sec.  Comparison  of  the  skin  friction  coefficients  are  shown  in  Figure 
2-3  for  freestream  velocities  of  32  and  242  ft/sec.  Even  though  the  values 
of  the  calculated  skin  friction  coefficient  overpredicts  the  experimental 
values,  the  equivalent  sand-grain  height  calculated  from  equation  (2-18)  is 
in  excellement  agreement  with  the  value  used  by  Healzer,  et  al .  These  values 
are  0.032  inches  for  the  predictions  as  compared  to  0.031  inches  suggested 
by  Healzer,  et  al .  based  on  Schl ichting' s  data.  Since  the  roughness  elements 
in  the  Healzer,  et  al .  experiments  were  touching  each  other,  it  was  necessary 
to  define  a  new  y  origin.  The  y  location  at  D(y)  =  0.99  D(o)  was  arbitrarily 
chosen  as  the  reference  level.  Velocity  profile  predictions  and  data  are 
compared  in  Figure  2-4.  Since  predicted  skin  friction  coefficients  were 
high,  predicted  velocity  profiles  were  normalized  by  the  experimental  edge 
values  of  u+.  Again,  it  appears  that  predicted  velocities  are  too  low  very 
near  the  wall.  Comparisons  to  the  rough-wall  predictions  for  momentum 
thickness  are  shown  in  Figure  2-5.  The  increase  in  the  values  of  0  due  to 
roughness  and  the  lesser  sensitivity  to  freestream  velocity  are  correctly 
predicted. 


No  rough-wall  skin  friction  data  seems  to  exist  for  3-D  flows. 
Therefore,  only  hypothetical  cases  could  be  considered.  Sample  cases  were 
run  for  a  Reynolds  number  representative  of  a  full-scale  ship  (RL  =  2  x  109 ) 
and  a  Reynolds  number  representative  of  a  wind  tunnel  model  (RL  =  lxlO7). 
The  sample  ship  was  the  same  elliptic  cross-section  body  a*-  the  sample  case 
of  K.  A  reference  smooth  wall  case  was  run  at  each  Reynolds  number  along 
with  roughnesses  defined  by  an  equivalent  sand-grain  height  or  roughness 
element  height,  base  diameter  and  center-to-center  spacing.  For  the 


roughness  calculations  it  is  also  necessary  to  input  a  characteristic  length 
(ship  length).  In  the  present  form  of  the  computer  code  only  spherical 
segment  roughness  elements  are  considered.  For  the  results  presented  here 
k  =  0.00083  ft.,  D(0)  =  0.00167  ft.,  and  i  =  0.004175  ft.  The  ship  lengths 
used  were  1000  feet  for  the  full  scale  model  and  10  feet  for  the  wind  tun¬ 
nel  model.  Plots  of  the  streamwise  momentum  thickness  and  the  crossflow 
angle  8  versus  arc  length  along  the  girth  at  X/L  =  0.84  are  given  in  Figures 
2-6  and  2-7  for  RL  =  1  x  107  and  in  Figures  2-8  and  2-9  for  RL  =  2  x  109. 

As  expected  the  momentum  thickness  increases  with  roughness.  However,  as 
shown  in  Figures  2-6  and  2-8,  this  increase  is  significantly  less  near  the 
keel.  This  phenomenon  is  noticeable  as  one  approaches  the  stern.  Roughness 
was  predicted  to  enhance  the  crossflow  angle  8  (see  Figures  2-7  and  2-9). 

This  is  especially  true  as  the  waterline  is  approached.  Intuitively,  this 
result  was  not  expected.  It  was  felt  that  surface  roughness  would  tend  to 
wash  out  three-dimensional  effects,  not  enhance  them.  Measurements,  or 
perhaps  more  exact  calculations,  are  needed  to  more  precisely  define  sur¬ 
face  roughness  effects  in  three-dimensional  flows. 

These  3-D  calculations  only  serve  the  purpose  of  indicating  that 
this  roughness  theory  has  been  implemented  and  debugged  in  the  computer  pro¬ 
grams  of  von  Kerczek  (1983)  corresponding  to  the  3-D  boundary  layer  theory  of 
K.  It  would  be  very  valuable  to  do  some  numerical  experiments  with  this  pro¬ 
gram.  One  set  of  tests  that  would  be  valuable  is  a  calculation  of  the 
drag  increment  ACD  obtained  by  calculating  the  skin  friction  drag  for  a 
smooth  and  a  rough  hull.  The  form  drag  of  Granville  (1974),  which  only 
depends  on  CD  and  the  potential  flow,  can  be  used  to  complete  the  determina¬ 
tion  of  ACg.  Although  no  detailed  3-D  roughness  boundary  layer  measurements 
seem  to  be  available,  there  probably  is  some  data  (possibly  for  full  scale 
ships)  for  which  ACD  has  been  determined.  For  example,  surely  trials  on  a 
new  ship  and  then  later  when  the  ship  has  aged  have  been  performed.  A  thor¬ 
ough  literature  search  for  such  cases  would  be  highly  desirable. 


b.  Experimental 


o  Experimental 
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Figure  2-3.  Rough  surface  friction  versus  x-Reynolds  number 
(experimental  data  from  Healzer,  et  al . ,  1974). 
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Figure  2-4.  Velocity  profiles  at  different  freestream  velocities 
(experimental  data  from  Healzer,  et  al . ,  1974). 
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Figure  2-5.  Rough  surface  momentum  thickness  versus  distance 
(experimental  data  from  Healzer,  et  al . ,  1974). 
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Figure  2-6.  Girthwise  distribution  of  streamwise  momentum 

thickness  8u  for  the  elliptic  cross-section  ship 
model  RL  =  1  x  10?. 
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Figure  2-7.  Girthwise  distribution  of  cross-flow  angle  B  for 
the  elliptic  cross-section  model,  RL  =  1  x  107- 
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Figure  2-8.  Girthwise  distribution  of  momentum  thickness  9ii 
for  the  elliptic  cross-section  ship,  RL  =  2  x  10- 
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Figure  2-9.  Girthwise  distribution  cross-flow  angle  3  fo 
the  elliptic  cross-section  ship,  RL  =  2  x  10 
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VISCOUS- INVISCID  INTERACTIONS 


3. 1  Introduction 


It  has  been  well  established  now  that  the  stern  boundary  layer  and 
near  wake  cannot  be  handled  by  classical  boundary  layer  theory  alone  (see 
Larson,  1981  and  Patel,  1983).  The  stern  boundary  layer  is  fairly  thick  and 
causes  a  considerable  disturbance  to  the  bare  body  potential  flow.  There 
are  two  currently  popular  ways  to  deal  with  the  stern  flow.  One  way  is  by 
using  boundary  layer  theory  with  viscous-inviscid  interactions  and  the  other 
way  is  by  abandoning  the  boundary  layer  equations  in  favor  of  the  "Parabol¬ 
ized  Navier-Stokes"  (PNS)  equations  (see,  for  example,  Hogan,  1983).  The 
latter  method  is  very  complex  and  relies  more  heavily  on  large  computers 
and  computing  budgets.  In  keeping  with  the  overall  simplicity  and  practi¬ 
cality  of  the  momentum  integral  method  it  is  desirable  to  develop  a  rela¬ 
tively  simple  way  to  do  viscous-inviscid  interaction  calculations.  It  is 
not  expected  that  such  a  method  will  be  as  accurate  as  PNS  calculations  (if 
they  ever  do  become  fully  developed)  or  even  viscous-inviscid  interaction 
methods  based  on  finite  differenced  full  higher  order  boundary  layer  equa¬ 
tions  (if  these  can  ever  be  made  to  work  for  ship  sterns),  but  at  least  a 
small  improvement  at  the  stern  of  the  basically  simple  and  practical  momen¬ 
tum  integral  method  seems  worthwhile. 

A  considerable  effort  has  been  expanded  to  develop  a  momentum 
integral-based  method  of  viscous-inviscid  interaction.  This  effort  was  at 
first  unsuccessful.  (For  completeness,  a  brief  outline  of  this  effort  is 
given  in  the  appendix.)  Unfortunately,  the  limited  budget  of  this  contract 
did  not  allow  completion  of  a  new  method  which  has  been  developed  and  is 
believed  to  hold  considerable  promise.  This  new  method  of  viscous-inviscid 
interaction  is  detailed  below.  The  method  is  based  on  the  boundary  layer 
momentum  integral  equations  and  slender  body  potential  flow  theory. 


A  word  should  be  said  here  in  defense  of  slender  body  potential 
flow  theory.  This  theory  has  often  been  cited  as  being  inadequate  for 
describing  the  pressure  distribution  on  bodies,  particularly  near  the  ends. 
This  is  only  partially  true.  However,  slender  body  theory  need  not  be  used 
to  calculate  the  boundary  layer  over  the  first  60  to  80  percent  of  the  body. 
There  the  boundary  layer  is  thin  and  does  not  siginficantly  interact  with 
the  body's  potential  flow.  Hence,  this  part  of  the  boundary  layer  can  be 
calculated  once  and  for  all  using  an  exact  potential  flow.  The  slender  body 
theory  only  comes  into  play  in  the  stern  region.  But  there  the  viscous-in- 
viscid  interaction  effectively  smooths  out  the  body,  i.e..  the  displacement 
body  is  much  smoother  and  more  gently  sloping  than  the  actual  body.  Hence 
there  is  no  reason  to  believe  that  slender  body  theory  cannot  accurately 
describe  the  potential  flow  over  the  displacement  body.  Furthermore, 
with  slender  body  theory  an  infinite  displacement  body  in  the  wake  can  be 
accounted  for  with  no  truncation. 

The  new  viscous-inviscid  interaction  theory  is  described  in  the 
next  two  subsections.  In  order  to  clearly  illustrate  this  theory,  the 
axisymmetric  case  is  treated  first.  Then  extension  to  three-dimensional 
ship  flows  is  outlined. 

3.2  Axisymmetric  Viscous-inviscid  Interactions 

The  starting  point  of  this  theory  is  the  following  physical  pic¬ 
ture  of  viscous-inviscid  interaction  employed  successfully  by  Huang,  et  al. 
(1976).  The  boundary  layer  develops  on  the  body  defined  by  the  radius  rQ(x), 
xe  [x  ,  xj,  where  x  is  the  axial  coordinate,  x  is  the  body's  nose  location 
and  xt  is  the  body's  tail  location.  The  boundary  layer  flow  is  driven  down¬ 
stream  on  the  body  by  the  inviscid  flow  corresponding  to  the  displacement 
body  defined  by  rQ(x)  +  6*(x)  cosy  where  6*(x)  is  the  effective  displacement 
thickness  and  S*(x)  0  for  x  -*■  «.  Huang,  et  al .  developed  a  method 

which  relies  on  calculating  the  boundary  layer  only  to  the  end  of  the  body 
(or  nearly  to  the  end,  say  x  =  b<xt).  The  extension  of  the  displacement 
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thickness  5*(x)  into  the  wake  is  handled  by  an  extrapolation  polynomial  which 
is  controlled  far  downstream  by  an  estimate  of  the  total  body  drag.  The 
iteration  procedure  to  solve  this  problem  is  developed  on  the  basis  that 
successive  calculated  displacement  bodies  have  the  same  pressure  distri¬ 
bution.  This  method  of  iterating  the  boundary  layer  and  potential  flow 
seems  to  work  fairly  well  even  though  it  does  not  really  converge  very 
rapidly.  The  method  seems  to  oscillate  about  a  mean  value  of  drag.  Huang, 
et  al . ,  simply  average  these  oscillations  to  obtain  a  fairly  accurate  final 
result.  The  method  outlined  below  makes  apparent  why  this  type  of  method 
behaves  in  the  way  Huang,  et  al.,  discovered.  It  is  noted  here  that  almost 
all  viscous- inviscid  interaction  methods  seem  to  have  such  convergence 
problems.  It  is  found  that  rather  drastic  under-relaxation  needs  to  be 
used  to  make  such  methods  converge  (see,  for  example.  Carter,  1978).  This 
convergence  problem  may  portend  serious  difficulties  for  3-D  flows.  Such 
flows  have  an  extra  degree  of  freedom  which  usually  results  in  more  diffi¬ 
cult  convergence  problems. 


The  starting  point  of  this  development  is  the  slender  body  theory 
of  Tuck  and  von  Kerczek  (1968).  In  polar  coordinates  with  r  denoting  the 
radial  and  x  the  axial  coordinates,  the  velocity  potential  is  given  explicit¬ 
ly  by  the  formula 
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4>  =  x  +  c(x)  Jin 


?0(x)r-  Js / dc(£) 


sgn(x-C)  £n  2 |x  -  £ 


(3-1) 


where  r  (x)  =  r  (x)  +  6*(x)  cosy  (r  (x)  =  r  (x)  in  the  case  of  no  viscous- 
o  o  0  0 

inviscid  interaction),  r  =  1  evaluates  $  on  the  body  rQ  and 


c(x) 


r0(x) 


d  y*) 

dx 


The  upper  limit  of  the  integral  is  xt  when  <$*(x)  =  0. 


(3-2) 


It  should  be  noted  that  if  the  body  has  a  blunt  tail  the  potential 
(3-3)  is  very  singular  there,  and  hence  the  boundary  layer  on  such  bodies 
separates.  This  could  prevent  any  viscous-inviscid  interaction  method  from 
converging  because  a  separated  flow  is  not  a  suitable  first  guess.  This 
can  be  overcome  by  artificially  extending  the  body  into  the  wake.  This  pro¬ 
cedure  can  be  carried  out  by  the  same  method  that  is  described  below  for 
extending  the  displacement  surface  into  the  wake.  Thus,  a  way  to  obtain  a 
suitable  starting  boundary  layer  needs  no  further  discussion. 

It  is  now  assumed  that  calculations  of  the  boundary  layer  have 
yielded  the  displacement  surface 


rQ(x)  =  rQ(x)  +  6*(x)  cosy  (3-3) 

up  to  x  =  b<xt  ( |  (x  -  b)  |  <<1 )  and  where  y  is  the  angle  between  the  tangent 
to  the  body  meridian  and  the  x-axis.  Then 


c(x) 


=  r  r 
o  o 


A)_ 

2tt 


(3-4) 


where  A  is  the  cross-sectional  area  of  the  displacement  body  and  primes 
denote  differentiation  with  respect  to  x. 


The  slender  body  theory  (3-1)  requires  that  c(x)  -*■  0  as  x  °°  in 
order  that  <J> <  °°.  Thus  it  is  best  to  work  directly  with  extrapolating  c(x) 
downstream  rather  than  rQ(x).  From  here  on  subscribts  b  denote  evaluation 
at  x  =  b. 


The  following  inverse  polynomial  extrapolation  has  suitable 
behavior  at  x  +  °°  and  can  be  matched  to  the  displacement  body  at  x  =  b: 

2tt  c'(x)  =  A"  ( x)  =  (A"  +  K)  (£^)2+3  -  K(^)3+3  (3-5) 

where  X  <  b,  6>0  and  K  are  adjustable  parameters.  Note  that  equation  (3-5) 
satisfies  the  two  requirements  that  A"(b)  =  A"^  and  A"(°°)  =  0.  By  integrating 
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equation  (3-5)  twice  and  matching  to  the  displacement  body  at  x  =  b,  the 
wake  displacement  body  r(x)  is  obtained.  Thus, 


and 
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(3-6) 
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where 


K '  -  01  f(b-x)  Ab +  (I+B> Abl 


(3-8) 


The  evaluation  of  equation  (3-7)  for  yields  the  relationship 
between  A^  and  A^,  namely 


a.  -  v  *  !(i-x)  Ai,  *  (b-x)2  Ab 


(3-9) 


This  extrapolated  wake  is  C2  continuous  (i.e..  A,  A',  A"  are  continuous)  at 
x  =  b.  Thus,  since  the  quantities  A^ ,  A^  and  A£  are  computed  by  boundary 
layer  theory  it  is  expected  that  wake  (3-7)  is  an  accurate  approximation  of 
the  wake  in  a  significant  neighborhood  of  the  point  x  =  b.  If  equation  (3-7) 
departs  significantly  from  the  asymptotic  shape  of  the  wake,  the  departure 
should  not  have  too  much  of  an  effect  on  the  potential  flow  velocities  up¬ 
stream  of  and  at  the  point  b. 


It  is  known  (see  Schlichting,  1979)  that  the  far  wake  displacement 
area  A^  is  related  to  the  drag  coefficient  of  the  body  by  the  equation 

*tcd-  (3-10) 

where  AQ  is  the  body  reference  area.  Hence,  equation  (3-9)  yields  the  rela¬ 
tionship 
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A0  CD  =  4Ab  +  f  (b‘A)  Ab  +  (1+3H'2+"B)  (b‘X)  Ab 


(3-11 


which  will  be  the  basis  of  the  iteration  procedure  for  calculating  the  vis- 
cous-inviscid  interaction. 

However,  before  describing  this  iteration  procedure  in  detail  it 
is  best  to  carry  equation  (3-11)  a  little  further.  The  axisymmetric  momen¬ 
tum  integral  boundary  layer  equations  are  the  following: 


f  +0(2*H)if  =  Hr0(x)  cf 


I?  *  o  Ij  fs  ({C0SY  +  ro)  F(G) 


( 3- 12a 


(3- 12b 


where  0  is  the  momentum  area,  H  =  A*/0,  A*  is  the  displacement  area  (note 
A*  f  A),  Q  =  GG,  G  is  Head's  shape  parameter,  s  is  arclength  along  a  body 
meridian,  cf  is  the  skin  friction  coefficient,  F  is  the  entrainment  rate  and 
U  is  the  potential  flow  speed  on  the  displacement  body.  See  K  for  the 
details  of  these  boundary  layer  equations.  The  arclength  s  and  the  axial 
coordinate  x  are  related  by 


ds  =  7  1  +  r' (x)2  dx. 


In  the  vicinity  of  the  tail  of  the  body 


V  0 


Cf  *  0 


ds  s  dx  (or  dx) 


H  -  constant  =  HL . 


(3-13 


Thus,  the  momentum  equation  (3- 12a)  reduces  there  to  the  equation 


A'  +  A(2+Hb)  U'/U  =  0. 


(3-14) 


By  evaluating  this  equation  at  x  =  b,  the  following  relationships  are 
obtained: 

(2+H.  )  U' 

A'  » - ^ ^  Ah 
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(3-15a) 
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( 3- 15  b ) 


Hence,  equation  (3-11)  can  be  rewritten  as 
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(b-X)2(2+Hb)[(3^Hb)  U’2  -  Ub  Ub 


( 1+3 ) ( 2+3 )  U2 


(3-16) 


Equation  (3-16)  clearly  brings  out  the  direct  and  seemingly  sensitive 
dependence  of  CD  on  Ub,  Ub  and  Uj\  In  fact,  it  seems  likely  that  a  good 
choice  of  A  is  such  that  (b-A)  =  h,  say,  and  that  probably  3  =  0(1).  Then 
the  third  term  in  equation  (3-16)  is  probably  negligible.  It  is  then  fairly 
easy  to  see  that  viscous-inviscid  iteration  schemes  will  oscillate  at  each 
iteration.  Equation  (3-16)  is  approximately  the  same  as 


(2+HJ 


floC0  =  4Ab<1-  (Sib 


up. 


(3-17) 


and  under  the  same  assumptions,  equation  (3-11)  is  approximated  by 


A  Cn  =  4  A.  +  |  A'  . 
o  D  b  6  b 


(3-18) 


3-7 


Assume  that  n  is  the  iteration  number.  If  at  iteration  n  the 
value  of  >  0,  then  it  will  produce  a  new  velocity  Ubn+1 >  0.  Then  if  this 
new  velocity  distribution  for  which  >  0  is  used  to  recalculate  the 

boundary  layer,  producing  A^j,  then  by  equations  (3-17)  and  (3-18)  (and 
also  simple  physical  reasoning)  it  is  evident  that  ^^<0.  Thus,  it  is 
apparent  that  viscous-inviscid  iteration  will  always  lead  to  oscillatory 
convergence  (or  divergence). 


It  is  now  worthwhile  to  state  the  iteration  algorithm  for  viscous- 
inviscid  interaction  calculations. 


Step  0:  rQ(x) 


Step  1: 


Step  2: 


p,x>  x"SxSx‘} 

L  o  x  >  x+  J 


If  body  is  very  blunt  at  the  tail,  possibly  causing  boundary  layer 
separation,  add  an  artificial  tail  to  the  body  using  wake  formulas 
(3-5)  -  (3-7)  with  b  =  x^  -  e,  e>0,  r  ^  =  rQ(xt  -  e)  and  c^  = 
c'(xt  -  e)  and  replacing  rQ(x)  =  0  with  rQ(x)  =  this  artificial 
tail  for  x  >xi.. 


Calculate  potential  flow  on  this  body  (either  exact  or  by  slender 
body  theory  or  a  combination  of  these). 


Calculate  boundary  layer  on  body  rQ(x)  to  x  =  b  to  get  A(x)  and 


A|j»  A^,  A^. 


Construct  wake  function  using  A^,  A^,  A£  and  calculate  C^  from 
formula  (3-11).  Is  CD  =  CD  from  previsou  iteration?  (if  on  first 
step  go  to  3)  Yes;  go  to  step  4.  No;  go  to  step  3. 


Step  3:  Calculate  U  by  slender  body  theory  potential  (3-1)  using  A(x)  on 
body /wake  from  steps  1-2.  Go  to  step  1. 


Step  4:  Calculate  CD  by  stress  integration.  If  Cn  is  same  as  wake  C 


D 
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problem  is  finished.  If  not,  adjust  X(or  0)  so  that  the  two  values 
of  Cn  are  equal  and  restart  at  step  2. 


Three-Dimensional  Viscous-Inviscid  Interactions 


The  beauty  of  the  theory  outlined  in  the  previous  subsection  is 
that  it  can  be  applied  with  only  minor  modifications  to  a  ship  hull.  The 
starting  point  of  this  development  is  the  representation  of  a  ship  hull  by 
conformal  mapping  functions  (von  Kerczek  and  Tuck,  1969)  and  the  general 
ship  slender  body  theory  (Tuck  and  von  Kerczek,  1968). 

The  ship  hull  is  given  by  r=l  in  the  equation 


■ro(x,0,r)  =  xi  +  y(x,B,r)j  +  z(x,6,r)k. 


In  equation  (3-19) 


N  ( 3-2n) 

y  +  iz  =  £  a  (x)c 


where 


and  an(x)  are  interpolation  functions  along  the  axis  of  the  ship. 

The  slender  body  velocity  potential  cj>(x,r,e)  is  obtained  directly 
in  terms  of  the  coefficient  functions  an(x)  by  the  following  formulas: 

(Tuck  and  von  Kerczek,  1968) 

N-l  c  (x)  P 

.  .  .  /..\  „  _  _  /  ^  n  -Zn 


<J>  =  X  +  CQ ( x )  £n  aj(x)c  -  H  -  C 

n=  1 

oo 

-  hf  d  cQ(C)  sgn(x-C)  An2|x-?| 


where  the  coefficients  cn(x),  n  =  0,  ....  N  are  obtained  from  the  formulas 


c  =  (y  +  Y  )  e 
n  'n  -n  n 


(3-22 


Y  =  E  (3-2ft)  a.(x)  a'  (x) 


(3-22 
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( 3- 22c ) 


The  potential  (3-22)  splits  into  two  important  and  distinct  parts. 
The  first  three  terms  of  equation  (3-22)  are  local.  They  only  depend  on 
the  particular  point  on  the  hull  at  which  the  velocity  is  evaluated.  It  is 
only  this  local  group  of  terms  that  contains  the  entire  non-axi symmetry  of 
the  hull  shape.  The  integral  in  equation  (3-22)  is  the  only  term  that  con¬ 
veys  the  global  aspects  of  the  hull  to  each  point  on  the  ship.  This  term 
is  exactly  the  same  as  in  the  axisymmetric  case.  The  function  cQ(x)  is  the 
x-derivative  of  the  cross-section  area  divided  by  2ir.  Thus,  the  effects  of 
the  wake  are  felt  on  the  hull  only  by  its  section  area  distribution,  not 
its  cross-section  shape.  This  feature  of  slender  body  theory  is  the  crucial 
element  that  makes  the  axisymmetric  interaction  theory  easily  extendable  to 
three  dimensions. 


After  the  boundary  layer  is  calculated,  a  displacement  thickness 
distribution  6*(x,e)  is  obtained.  Since  this  displacement  thickness  is 
supposed  to  represent  a  distance  normal  to  the  hull,  the  equivalent  distance 
normal  to  the  cross  section  out  to  the  effective  displacement  surface  is 
given  by  6*c(x,6)  where 

6*(x,9)  =  6*(x  ,0) (TT  •  n) ,  (3-23) 

V*  V# 

where  n  is  the  unit  normal  vector  to  the  hull  and  nc  is  the  unit  normal 
vector  to  the  cross-section.  Furthermore,  near  the  stern  the  value  of  S*c 
may  be  so  large  that  if  this  distance  is  plotted  normal  to  the  cross-section, 
the  neighboring  displacement  distances  would  cross,  as  sketched  in  the  figure 
below. 
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Schematic  diagram  showing  inversion  of  outer  edge  of  boundary  layer 
when  plotted  as  a  function  of  distance  along  the  normals  (points 
denoted  by  dots,  •)  as  opposed  to  the  distance  along  the  coordinates 
r  (points  denoted  by  x's). 

This  difficulty  is  easily  overcome  by  plotting  the  displacement  thickness 
along  the  r-coordinate  curves  (8  =  constant)  in  the  cross-sections.  The 
cross-section  equation  (3-21)  is  a  conformal  mapping,  hence  the  0  =  constant 
curves  are  locally  orthogonal  to  the  cross-sections  (r  =  constant  =  1). 

Thus,  if  hr  denotes  the  metric  coefficient  of  the  r-coordinate,  then  at 
each  station  x  the  displacement  surface  "rjj(x,0)  is  obtained  by 


where 


"rd(x,0)  =  rQ(x,0,r6) 
r6(x,0) 
<5*c(x,0)  ~J 


hr(x,0)  dr. 


(3-24) 


(3-25) 


For  the  case  where  <5*c  is  small  enough  not  to  cause  the  problems  exhibited 
in  the  figure,  the  displacement  surface  can  be  obtained  simply  as 


"^(x,©)  =T(x,0,l)  +  6*c(x,0)"nc. 


(3-26) 
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Once  the  surface  "r^(x,0)  for  xc  [a.b]  has  been  constructed,  then  the 
cross-sections  of  this  surface  are  mapped  onto  the  unit  circle  (using  the 
von  Kerczek-Tuck,  1969,  algorithm)  and  a  new  set  of  gapping  coefficients 
an(x)  is  obtained.  From  these  coefficients  the  new  set  of  coefficients 
cn(x)  are  computed  by  formulas  (3-22). 

Now  the  section  area  derivative  cQ(x)  must  be  extrapolated  to 
infinity  to  complete  the  data  for  the  potential  (3-21).  Again  formulas 
(3-5)  -  (3-11)  are  used  for  the  wake.  The  same  algorithm  that  was  used  for 
the  axisymmetric  wake  can  now  be  used.  Of  course,  the  fully  3-D  boundary 
layer  calculation  method  of  K  (or  any  other  3-D  boundary  layer  method)  must 
be  used  to  obtain  the  new  effective  displacement  thickness  6*(x,0). 

The  beauty  of  this  method  is  that  the  actual  cross-section  shape 
of  the  wake  need  not  be  known  as  long  as  the  potential  flow  velcoity  on 
the  wake  displacement  surface  does  not  have  to  be  computed.  This  is  pre¬ 
cisely  what  is  achieved  by  extrapolating  the  wake  cross-section  area. 

Slender  body  theory  makes  this  possible.  It  should  be  emphasized  that  this 
theory  does  not  assume  that  the  wake  is  axisymmetric.  Use  of  slender  body 
theory  for  the  potential  flow  simply  obviates  the  necessity  to  know  the 
cross-section  shape  of  the  wake.  This  method  for  doing  three-dimensional 
viscous-inviscid  interaction  seems  extremely  promising.  Implementation  of 
the  method  to  bring  it  to  fruition  is  planned. 


COMPUTATIONAL  IMPROVEMENTS  TO  THE  BASIC  MOMENTUM  INTEGRAL  METHOD 
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It  was  found,  in  the  course  of  several  applications  of  the  method 
K,  that  the  integration  of  the  boundary  layer  equations  broke  down  prematurely 
near  the  stern.  This  breakdown  was  traced  to  .wo  causes.  The  first  cause 
was  an  inadequate  handling  of  the  equations  on  the  symmetry  plane.  The 
second  cause  was  the  requirement  in  the  original  program  that  the  boundary 
layer  equations  be  integrated  with  an  axial  step  size  conforming  to  the 
potential  flow  panels.  The  remedies  to  these  di fficiencies  which  required 
considerable  effort  are  described  below. 

The  original  computer  program  of  K  (see  von  Kerczek,  1983)  imple¬ 
mented  the  symmetry  plane  conditions  by  appropriate  parity  of  the  finite 
difference  formulas  on  these  planes  (the  load  waterline  and  keel).  It  was 
found  that  this  kind  of  implementation  of  symmetry  plane  conditions  allows 
the  boundary  layer  solution  off  the  symmetry  plane  to  affect  the  sym¬ 
metry  plane  solution.  But  boundary  layer  theory  indicates  that  the  flow 
along  a  symmetry  plane  is  completely  independent  of  the  flow  elsewhere  on 
the  body  (unless  viscous-inviscid  interactions  are  present).  Hence,  the 
symmetry  plane  boundary  layer  equations  must  be  solved  independently  of  the 
rest  of  the  boundary  layer  and  the  results  must  be  used  as  lateral  boundary 
conditions  to  the  full  three-dimensional  equations.  This  method  of  handling 
the  symmetry  planes  has  now  been  implemented  in  the  method  of  K  and  results 
in  improved  predictions  of  the  boundary  layer. 

The  original  computer  program  for  the  method  of  K  used  axial  step 
sizes  fixed  to  the  potential  flow  panels.  The  potential  flow  panels  cannot 
be  made  small  enough  at  a  ship  stern  to  accommodate  the  requirements  of  the 
boundary  layer  equations  unless  an  enormously  expensive  potential  flow  cal¬ 
culation  can  be  tolerated.  This  problem  was  alleviated  by  writing  a  pre¬ 
processor  that  interpolates  the  potential  flow  data  onto  an  arbitrary  grid. 
This  decouples  the  numerical  algorithm  for  solving  the  boundary  layer  equa¬ 
tions  from  the  potential  flow  panels.  The  boundary  layer  equations  in  K  can 


now  be  integrated  on  very  fine  step  sizes  near  the  stern  of  a  ship.  This 
greatly  improves  the  performance  of  the  code  of  K.  For  example,  the  bound¬ 
ary  layer  can  now  be  computed  to  within  2%  of  ships'  length  to  the  stern  of 
the  semi-elliptic  hull  of  Huang,  et  al.  (1976).  Comparison  of  such  calcu¬ 
lations  with  the  experimental  data  of  Huang,  et  al.  is  not  shown  here 
because  the  experiments  had  insufficient  detail  to  allow  the  discernment  of 
an  appropriate  set  of  initial  conditions  to  start  the  calculations.  The 
calculations  that  were  made  used  rough  guesses  for  initial  conditions,  and 
thus  the  results,  though  qualitatively  correct,  are  not  comparable  to  the 
experimental  data. 

It  is  planned  that  the  SSPA  model  boundary  layer  (see  K)  will  be 
recalculated  and  compared  to  experiment  with  the  improved  code  described 


APPENDIX:  UNSUCCESSFUL  METHOD  FOR  CALCULATING  VISCOUS-INVISCID 

INTERACTING  AT  A  STERN 

A  brief  outline  is  given  here  describing  an  unsuccessful  idea  for 
calculating  the  stern  flow  of  ships.  This  outline  is  given  mainly  to  illus¬ 
trate  the  difficulties  that  are  inherent  in  trying  to  join  computationally 
the  boundary  layer  and  near  wake  calculations.  The  difficulties  lie  mainly 
in  trying  to  match  up  two  geometrically  dissimilar  regions.  It  was  shown 
in  Section  3  that  a  viscous-inviscid  interaction  scheme  will  be  very  sensi¬ 
tive  to  the  potential  flow  velocity  and  its  first  two  derivatives.  This 
forebodes  difficulties  for  any  method  that  attempts  to  actually  match  up 
surfaces  in  three  dimensions. 

The  idea  in  the  first  attempt  to  do  viscous-inviscid  interactions 
was  to  model  the  displacement  effects  by  a  distribution  of  extra  sources  on 
the  body  and  its  extension  in  the  wake.  This  attempt  was  carried  out  on  the 
Groves,  et  al .  (1982)  semi-el  1 iptic  body. 

The  body  was  truncated  at  a  station  xt  -  c,  0<e«  1,  and  an  ellip¬ 
tic  cross-section  tail  cone  was  attached  there.  This  tail  cone  extended  a 
distance  downstream  determined  at  x  =  x^  -  e  by  the  slopes  of  the  waterline 
and  keel  of  the  actual  body. 

The  first  iteration  calculated  the  full  three-dimensional  bare 
body  boundary  layer  up  to  station  xt  -  e.  At  station  xt  -  e  onwards  the 
wake  was  computed  on  the  waterline  and  keel  of  the  extended  tail  cone. 

The  centerline  pressure  distribution  behind  the  body  without  the  tail  cone 
was  used  for  the  wake  calculation.  The  tail  cone  only  serves  the  purpose 
of  providing  a  continuation  of  the  coordinate  system  off  the  end  of  the  body. 

This  boundary  layer  and  wake  calculation  then  provides  the  displace¬ 
ment  surface  <5*(x,0)  (using  the  coordinates  of  section  3.2  say).  Extra 


A-l 


sources , 


o  =  U  ( A- 1 ) 

where  s  is  the  arclength  along  the  streamlines,  were  then  distributed  on 
the  body  and  its  extension.  The  increment  to  the  original  exact  bare  body 
potential  flow  was  then  computed  by  slender  body  theory.  Using  this  new 
potential  flow,  the  boundary  layer  and  wake  were  recalculated. 

The  method  failed  because  of  extremely  erratic  behavior  of  the 
potential  flow  and  the  extra  source  strengths  a  in  the  neighborhood  of  the 
tail  of  the  body.  Examination  of  the  details  of  this  method  did  not  lead 
to  optimism  towards  a  continuation  of  this  technique,  so  it  was  abandoned. 

The  lesson  learned  here  is  that  it  is  very  tricky  to  actually  con¬ 
tinue  a  three-dimensional  body,  even  a  displacement  body,  into  the  wake. 

The  reason  is  that  the  potential  flow  is  so  sensitive  to  inflection  points 
(i.e.,  in  naval  architectural  terms,  the  fairness  of  the  body)  that  it  is 
too  difficult  to  extend  surfaces  in  three  dimensions.  The  slender  body 
method  outlined  in  section  3.2  only  requires  that  a  curve  be  continued 
downstream,  namely  the  section  area  curve.  This  is  the  same  as  what  is 
required  in  the  axisymmetric  case,  and  this  has  been  shown  to  work  by 
Huang ,  et  al .  (1976) . 
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